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Abstract 

In general a decay with a missing (not detected) particle can not be fully reconstructed apart from a few 
exceptions. For example, if the momentum of the decaying particle is known or if the missing energy in an event is 
measured precisely, then the missing particle 4-momentum can be determined. Here a new method is proposed that 
utilizes additional information about the topology of a decay: The direction from the primary to the secondary 
vertex combined with momentum conservation allows the determination of the missing particle momentum up 
to a twofold ambiguity. The semileptonic decay of the _B° mesons is considered as an example to illustrate this 
method and to compare its performance against conventional approaches. 
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1. Introduction 

Generally it is considered that the decay of a particle, where one of the final state products is not 
detected, can not be fully reconstructed. There are a few well known exceptions. If the momentum of 
the decaying particle is known and all but one of the decay products are detected and reconstructed 
then, obviously, it is possible to determine the missing particle's 4-momentum. Another possibility to 
measure the momentum of missing particle is based on the hermeticity of the detector where the missing 
energy is measured sufficiently well. An approximate method is used in studies of B®B® oscillations 
with semileptonic decays where the neutrino is a missing particle. The reconstruction of the B® meson 
momentum is a key ingredient in these studies. For this purpose, a fc-factor is introduced which is the ratio 
between the momentum of the measured decay products and the original B® momentum. The /c-factor 
is obtained from Monte Carlo (MC) studies. The use of the fc-factor induces a large error on the proper 
time of the B®. 
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In this paper we point out that by including additional topological information the neutrino 4- momentum 
can be reconstructed up to a twofold ambiguity. The study of B®B®- oscillations in the semileptonic de- 
cay of the B° meson at hadron colliders is used to illustrate the method of the neutrino 4-momentum 
reconstruction. 



2. Kinematics of the Semileptonic B° Decay 

Consider the decay B® — > Dj£ + v where the Dj meson decays into a final state with charged hadrons 
and all hadrons are reconstructed. The kinematics of such a decay can be described by four equations in 
three dimensions: 

y/m% + P% = ^m\ Dsi) + P\ DJ) + |P„| (1) 

Pb=P(d.i)+P* (2) 

In this system of equations there are six unknown variables: Pg, P*, where i stands for x, y and z. 
Since there are only four equations, the system is undetermined. 

The flight direction of the B® meson, V b , obtained from the primary and secondary vertex, provides 
the necessary information to solve for the neutrino momentum. Without loss of generality one can work 
in a two dimensional coordinate system: one axis (||) is defined by the B® flight direction V b, the other 
one (_L) is defined in the plane spanned by the two vectors Pd s i and V b and is orthogonal to the B® 
flight direction. 9 is the angle between the flight directions of the B a s meson and the D s £ system. Then 
one gets: 
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The momentum components of the D s l system and of the neutrino orthogonal to the B® flight direction 
are of equal magnitude and opposite sign: 

P { D 3 e) = -P v ± (9) 

Therefore the neutrino momentum component along the B® flight direction can be obtained with a 
straightforward calculation: 

Pj = -a ± (10) 
where 
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* U-^D^y {D s l)> \ r (D s l)> (D.l) 

Eq. 10 leads to two solutions, only one of them is correct. In an experimental environment there might 
be no solution of this equation if the radicand r is negative due to finite vertex and momentum resolutions. 
In the following we illustrate that despite these facts the method still provides competitive results with 
respect to conventional approaches. 



3. MC simulations 



To verify that the proposed method works not only in ideal conditions without experimental errors, a 
MC simulation has been developed to study B®B® - mixing in the semileptonic decay mode. The PYTHIA 
V6.227 [1] package has been used as an event generator with a center-of-mass energy of 14 TeV. We have 
generated 45000 events each containing two b quarks. One of them hadronizes to a B° meson and decays 
according to 

B° s -£>,p+i/ M) D s ^cPtt-, <j>^K+K-; (13) 
the second b quark hadronizes according to PYTHIA into a &-hadron and decays unconstrained. The 
decay in (13) is called 'signal' decay to distinguish it from the other 6-hadron decay. 

All hadrons from the signal decay are required to have pr > 1 GeV/ c, the muon is required to have 
pr > 3 GeV/ c. The track parameters and vertex positions (primary and secondary) have been smeared 
according to Gaussian distributions with the following parameters. The momentum uncertainty has been 
simulated by smearing the pscudorapidity with a v — 5.8 x 10~ 4 , the angle <j> with = 0.58 mrad and 
the inverse transverse momentum with <7(i/ PT ) = 0.013(Ge^/c) _1 . The primary vertex has been smeared 
with u xy = 20 fxm in both, x- and y-direction, the secondary vertex has been smeared with <j\\ =70 /im 
in flight direction of the B® and a± = 10 pm in the perpendicular direction. 



4. Proper time reconstruction 

The most important ingredient in the measurement of the B® oscillation frequency is the proper time 
calculated in the transverse plane 1 as follows: 

L xy m(B° s ) 

Pt(B°) ' (i4j 

where L xy is the flight path length, m(B®) and pr(B®) are the mass and transverse momentum of the 
B® meson respectively. 

In the semileptonic B® meson decay, the neutrino is not detected and the pt of the B® meson cannot 
be measured directly. Therefore a correction factor, derived from MC simulations, is introduced to scale 
the measured transverse momentum of the D s l system. Eq. 14 is modified as follows 

L xym (B°) 

CT ~ Pt(D s £) Xk > (15) 



1 Here the transverse plane is perpendicular to the beam direction, as is usual in collider experiments. 
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with the fc-factor estimated from MC simulations and calculated as 

Pt{D s £) 



k = 



Pt(B°) 



(16) 
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Fig. 1. Probability density for the fc-factor. 

This fc-factor (see Fig. 1) introduces a significant error on the B® momentum and, hence, on the proper 
time as well. To reduce the error by an average fc-factor, it is calculated in bins of m^ B iy We illustrate 
in Fig. 2 the B® meson transverse momentum resolution as a function of the invariant mass of the 
D s £ system for the fc-factor and the neutrino reconstruction method. It is evident that the resolution 
obtained with the neutrino reconstruction method is substantially better than with the fc-factor method. 
Note that the neutrino reconstruction method is sensitive to a± , therefore two different performances for 
<T_l_ = 10 /im and a± = 30 /im are shown. 
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Fig. 2. B® momentum resolution calculated with fc-factor and neutrino reconstruction methods. 

Fig. 3 illustrates the proper time resolution obtained with the fc-factor method (a) and the neutrino 
reconstruction method (b). For the latter method the closest to the true value solutions are filled in the 
histogram. The distributions are fitted with two Gaussian, the average width a is determined according 
to 



n w n 1 w w 



(17) 
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where a n (cr w ) and N n (N w ) are the width and normalization of the narrow (wide) Gaussian, respec- 
tively. For the fc-factor method a n — 100 fs (a w — 338 fs) and N n = 2700 (N w = 730) that gives an 
average a — 132 fs. For the neutrino reconstruction method a n = 77 fs (a w — 193 fs) and N n = 2300 
(N w — 660) that gives an average a = 91 fs. The neutrino reconstruction method provides a proper time 
resolution which is better than in the reconstruction done with the fc-factor method. 
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Fig. 3. Proper time resolution, (a) fc-factor method, (b) neutrino reconstruction method. 



5. Amplitude Fit 

A standard way to search for oscillations is the amplitude method described in Ref. [2]. Briefly it 
consists of the following steps. Candidates for the decay (13) are split into two samples. Events with the 
same flavor at production and decay time comprise the unmixed sample, events with the opposite flavor 
constitute the mixed sample. The two cases are distinguished by tagging the flavor of the B® meson at 
the production and the decay. Tagging at the production can be achieved by identifying the flavor of the 
second &-hadron in the event; tagging at the decay can be accomplished by identifying the charge of the 
Bg meson decay products. 

In each event, the proper decay time of the B® mesons is reconstructed and the two event samples are 
used to define the time-dependent asymmetry 

P — P 

z,\ 1 unmix 1 mix a / a ,\ 

a(t) = — (x Ax D x cos(Am s i), (18) 

* unmix i ^mix 

where D is a global dilution factor accounting for background, miss-tagging and proper-time resolution 
and A is the amplitude. P unm ix and P m ix are the time-dependent probability distribution functions for 
mixed and unmixed B® meson decays, respectively. 

In the fit to the asymmetry distribution, the oscillation frequency Am s is fixed, leaving the amplitude 
A as a free parameter. A scan over Am s is performed starting from zero. If the fixed value of Am s is 
consistent with the true one, the fitted amplitude will be equal to unity, else it is consistent with zero. 
The error on the amplitude value A is calculated according to the formula 

i r^r~ s + b Am >t 

aA = T^2w x \lsTB: x — xe 2 ' (l9) 
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where W is the mistagging probability, S the number of signal, B the number of background events, and 
at the proper time resolution. The amplitude method is sensitive to the tested value of Am s if a a is small 
compared to unity. The value of Am s at which 1 = 1.645 x a a is quoted as sensitivity limit. 

In the following we assume a mistagging probability of 40%, a number of 45000 signal events for 
semileptonic B® decays and a signal to background ratio of 1:1. The simulated oscillation frequency is 
taken as Am s = 17.25 ps _1 [3,4]. 

The result of the amplitude fit for the fc-factor method with the sensitivity curve is presented in Fig. 4a. 
The sensitivity of the method for the processed number of signal events is about 17 ps _1 . 

On Fig. 4b one can see the result of the amplitude fit done on the events reconstructed with the neutrino 
reconstruction method. The sensitivity to Am s is about 21 ps" 1 . In this case there is a clear peak in the 
amplitude fit at the input value of Am s — 17.25 ps _1 . 




6. Discussion 



On the one hand the use of the neutrino reconstruction method improves the momentum resolution of 
the B® meson. On the other hand it degrades the signal to background ratio and reduces the signal sample. 
First, since solving the system of equations, we end up with a quadratic equation with two solutions. One 
solution is the 'true', the other is the 'wrong' one. In such a way the number of background events is 
increased by the number of wrong solutions from the signal as well as from the background sample. 
Second, again due to the quadratic equation and the realistic resolutions, a certain fraction of events has 
a negative radicand r and, hence, there is no solution for the B® meson momentum. These events are 
excluded from the analysis (both from signal and background samples). 

With the resolutions used in this analysis the number of signal events is decreased by a factor of 2 
roughly (the same is true also for background) and the number of background events is increased by 
a factor of 3 (due to the 'wrong' signal and background events). Hence, in the neutrino reconstruction 
method the signal to background ratio is not 1:1 (like in the fc-factor method) but 1:3. Despite these facts 
the sensitivity of the neutrino reconstruction method is still at higher values of Am s compared to the 
fc-factor method. 
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Also the quality of both methods with respect to the primary (Fig. 5a) and secondary (Fig. 5b) 
vertex resolution is investigated, when all other resolutions are kept constant. From these figures one 
can conclude that the neutrino reconstruction method is more powerful than the fc-factor method if 
<t_i_ < 45/im and a XtV < 45/xm. For a typical hadron collider detector resolutions better than quoted above 
are achievable [5,6,7,8]. 
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Fig. 5. Sensitivity as a function of (a) the primary vertex resolution cr x ,y and (b) the secondary vertex resolution aj_ for 
neutrino reconstruction and fc-factor methods. 



7. Conclusion 



In this paper we have shown that the full neutrino reconstruction is possible in decays like the semilep- 
tonic B° s decay mode where all but one final state particle arc measured. For this one has to use additional 
topological information, in our example the direction of the B® meson momentum. The search for the 
B®B° S oscillations has been used as an illustration of the method. For the verification of our procedure 
we have used typical resolutions of hadron collider detectors. The sensitivity obtained with the neutrino 
reconstruction method is much higher than with the convenient method used for this decay mode. 

Finally, we would like to stress that the proposed method can be used not only in semileptonic B decays 
but in some other cases where the known topology of a decay can compensate for the incompleteness of 
kincmatical information. 
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